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Abstract 

We study complex potentials and related non-diagonalizable Hamiltonians with special 
emphasis on formal definitions of associated functions and Jordan cells. The non-linear 
SUSY for complex potentials is considered and the theorems characterizing its structure 
are presented. We define the class of complex potentials invariant under SUSY transfor- 
mations for (non-)diagonalizable Hamiltonians and formulate several results concerning 
the properties of associated functions . We comment on the applicability of these re- 
sults for softly non-Hermitian PT-symmetric Hamiltonians. The role of SUSY (Darboux) 
transformations in increasing/decreasing of Jordan cells in SUSY partner Hamiltonians 
is thoroughly analyzed and summarized in the Index Theorem. The properties of non- 
diagonalizable Hamiltonians as well as the Index Theorem are illustrated in the solvable 
examples of non-Hermitian reflectionless Hamiltonians . The rigorous proofs are relegated 
to the Part II of this paper. At last, some peculiarities in resolution of identity for discrete 
and continuous spectra with a zero-energy bound state at threshold are discussed. 

1 Introduction 

Quantum Physics of open systems often deals with incomplete information on the influence 
of an environment and can be adequately described by non-Hermitian Hamiltonians with 
a non-positive imaginary part. This kind of effective description has been employed in 
Condensed Matter, Quantum Optics and Hadronic and Nuclear Physics pQ - [1] for many 
years. Non-self-adjoint operators were also under mathematical investigations [5j EJ [7] 
and recently interesting examples of non-Hermitian effective Hamiltonian operators have 
been found for the quantum many-body equations [8] . 

The PT-symmetric Quantum Mechanics proposed in [91 [10] and developed in |10|-|14j 
and its pseudo-Hermitian generalization |15} HB[ [T7] describes a variety of non-Hermitian 
Hamiltonians with real spectrum (but not all Hamiltonians with real spectrum are PT- 
symmetric [18\ I19j). There is a progress in understanding some non-Hermitian but PT- 
symmetric Hamiltonians in terms of Krein spaces [20] . This kind of Quantum Mechanics 
has attracted much interest as it may open the way to give a solid probabilistic inter- 
pretation of non-Hermitian dynamics by means of a positive pseudo-norm |11} 112] . PT- 
symmetry endows with a physical meaning the energy spectrum of some Hamiltonians 
formally unbounded from below (9j [10]. The latter possibility for anharmonic oscillators 



with potentials unbounded from below was observed long ago [21] [22] but only recently 
has been associated with a PT-symmetry |24| . 

For complex, non-Hermitian potentials the natural spectral decomposition involves 
biorthogonal states [25j. Moreover the Hamiltonians may not be diagonalizable [26J but 
can be reduced only to a quasi-diagonal form with a number of Jordan cells [15| . This fea- 
ture appears at level crossing which, in fact, occurs under specific circumstances in atomic 
and molecular spectra J2f5] and Optics |27[ [25] as well as in PT-symmetric quantum sys- 
tems |23[ [29] 130] , There are also certain links [31J to the occurrence of non-Hermitian 
degeneracies for essentially Hermitian Hamiltonians where the description has been de- 
veloped for complex eigenvalue Gamow states (resonances) unbounded in their asymp- 
totics and, in general, not belonging to the Hilbert space of physical wave functions. On 
the contrary, in what follows we examine non-Hermitian Hamiltonians with normaliz- 
able bound and associated states. The subtleties of biorthogonality (the phenomenon of 
" self-orthogonality" [28] ) in resolution of identity and in definition of quantum averages of 
observables have been thoroughly analyzed in our paper [32J . 

We find it certainly interesting and important to investigate the possible ways for 
quantum design of such non-Hermitian quantum systems and in particular to extend 
the methods of non-linear SUSY algebra [33]-[6l] in order to keep under control the 
emerging of non-diagonal parts of those systems. In making a link to PT-symmetric 
systems we restrict ourselves with a soft type of non-Hermiticity when the real part of 
a potential dominates over the imaginary one at both infinities and asymptotically such 
a potential remains bounded from below. Respectively the energy spectrum of a related 
system contains a number of bound states and possibly a continuum part bounded from 
below. Thus having in mind the SUSY quantum design one can, for instance, think of a 
chain of complex Hamiltonians produced by Darboux-Crum transformations from a real 
one as a good representative of the class of softly non-Hermitian systems. The general 
relations and theorems presented in Sections 2-4 certainly hold also for PT-symmetric 
potentials with fixed asymptotics of ratio of imaginary and real part (semihard non- 
Hermiticity), say, for potentials with leading asymptotics Xx 2n (ix) e ; A > at infinities 
provided that the boundary conditions for eigenvalue problem do not require to move 
to complex coordinates [101 23J (z. e. for |e| < 1) . However the Lemmas and Theorem 
of Sections 5 and 6 are proven for softly non-Hermitian potentials and we pay hopes to 
extend them also on semihard non-Hermitian potentials in a nearest future. 

We start in Sec. 2 with the definitions and a summary of properties of non-Hermitian 
diagonalizable Hamiltonians and introduce the relevant biorthogonal expansions. Then 
we consider non-diagonalizable non-Hermitian Hamiltonians with discrete spectrum and 
finite-size Jordan cells and discuss the choice of the biorthogonal basis with diagonal 
resolution of identity. The novel result of this section is the proof that a biorthogonal 
basis always exists which is made of a set of eigenfunctions and associated functions of the 
initial Hamiltonian and a set of their complex conjugates for the Hermitian conjugated 
Hamiltonian. Moreover if the Hamiltonian is PT-symmetric and this symmetry is not 
spontaneously broken on states (eigenvalues are real) then the elements of direct and 
conjugated bases are related by PT-reflection. In Sec. 3 the origin of non-diagonalizable 
Hamiltonians is clarified to be level confluence. 

Non-linear SUSY in QM is summarized and extended to complex potentials in Sec. 4 
with an emphasis to the possibility of conservation of PT-symmetry. Herein the impor- 
tant theorem on the polynomial structure of SUSY algebra with transposition symmetry 
as well as the strip-off theorem describing the minimization of the differential order of 
intertwining operators are adapted to the complex potentials. These theorems involve 
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the zero-mode subspaces of supercharge components - intertwining operators and their 
mapping by Hamiltonians - matrices S. They are well compatible with PT-symmetry 
(if any). The relationship between superpotentials and Wronskians involving associated 
functions is discussed. 

In Sec. 5 we present the class of complex potentials invariant under SUSY transfor- 
mations for (non-)diagonalizable Hamiltonians: this class covers the systems with soft 
breaking of Hermiticity and essentially real continuum spectrum (if any). For this case 
we formulate several results characterizing the normalizability of associated functions at 
+00 and/or — 00. The necessary conditions for SUSY transformation functions are found 
to provide a pre-planned Jordan structure of a SUSY partner Hamiltonian. These re- 
sults allow to unravel the relation between Jordan cells in SUSY partner Hamiltonians 
the latter being described by the Index Theorem in Sec 6. It represents the main re- 
sult of the present paper. The Index Theorem relates the dimensions of Jordan cells of 
super-partner Hamiltonians at any energy level with characteristics of intertwining oper- 
ator kernels (matrices 5)and in fact exhaustively describes the quantum design options 
for softly non-Hermitian Hamiltonians. Needless to say that the latter theorem is also 
compatible with PT-symmetry when non-Hermiticity , for instance, is introduced into 
P-even potentials by shifting of coordinates into complex plane (see examples in [29]) . 
The illustration of properties of non-diagonalizable Hamiltonians as well as of the Index 
Theorem is thoroughly performed in Sec. 7 by the solvable example of non-Hermitian 
reflectionless Hamiltonians originated by SUSY transformations from the free particle 
Hamiltonian. The arising of non-diagonalizability is illuminated by an exactly solvable 
system with two coalescing bound states. In Conclusions we outline possible peculiari- 
ties of non-Hermitian Hamiltonians with continuous spectrum. The approaching to the 
continuum threshold yields more subtle problems with normalizable eigen- and associated 
functions in continuum which may have zero binorm. As a consequence it may cause 
serious problems with the resolution of identity investigated in detail elsewhere [32j . 

All the new results on the structure of non-diagonalizable SUSY Hamiltonians pre- 
sented in this (part of) paper are rigorously proved and justified in the accompanying 
(second part of) paper [65] . 

2 Non-Hermitian diagonalizable vs. non-diago- 
nalizable Hamiltonians and biorthogonal expansions 

In our paper we deal with complex one-dimensional potentials V(x) 7^ V*(x) and re- 
spectively with non-Hermitian Hamiltonians h of Schrodinger typa^J, defined on the real 
axis, 

h = -d 2 + v{x), (i) 

which are assumed to be symmetric or self-transposed under the * - transposition op- 
eration, h = h . The notation d = d/dx is employed. Only scalar local potentials will 
be analyzed which are obviously symmetric under transposition (for some matrix non- 
diagonalizable problems, see [661 EZj). Taking into account possible applications in PT- 
symmetric QM we specify complex potentials to give a semihard non-Hermiticity when 
the real part ReU is bounded from below and the ratio ImU/ReU remains finite and suf- 
ficiently small for large x — > ±00 . In this case the eigenvalue problem can be safely posed 

1 Conventionally the system of units m = 1/2, H = c = 1 will be used with dimensionless energies, momenta 
and coordinates. 
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keeping the boundary conditions at x ± oo on the real axis. Later on, in last two sec- 
tions we restrict ourselves with softly non-Hermitian potentials with vanishing asymptotic 
ratios ImF/ReV = o(l) . 

Let us first define a class of one-dimensional non-Hermitian diagonalizable Hamiltoni- 
ans h with discrete spectrum such that: 

a) a biorthogonal system {\ijjn}, IVVi)} exists, 

h\lp n ) = X n \lpn)i h^\tjj n ) = \* n \lpn), {i>n\lpm) = (i>m\i>n) = <W, (2) 

b) the complete resolution of identity in terms of these bases and the spectral decom- 
position of the Hamiltonian hold (in the case of PT-symmetric potentials the necessary 
conditions for that are formulated in |20j), 

I = ^2\^n}(lpn\, h = } ] X n \lp n )(lp n \. (3) 
n n 

In the coordinate representation, 

1p n (x) = {x\lp n ), Tpn(x) = (4) 

the resolution of identity has the form, 

5(x - x') = (x'\x) = ^ Mx'Wn(x). (5) 

n 

The differential equations, 

hipn = Ktpn, h^ n = A* 4> n , (6) 

and the fact that there is only one normalizable eigenfunction of h for the eigenvalue A n 
(up to a constant factor), allow one to conclude that 

4>n(x) = a n ip n (x), a n = Const / 0. (7) 

Hence the system {\i/j n ), IVVi)} can be redefined 

Hn) - -L^n), \4>n) -» V<$n), (8) 



so that 

4>*(x) = ip n {x), I ip n (x)ip m (x) dx = 5 nm . (9) 



We stress that the non- vanishing binorms in Eq.([9]) support the completeness of this basis, 
i.e. the resolution of identity, 

5{x - x) = S ^1p n {x)^ n (x). (10) 

n 

Indeed if some of the states in Eq. (fTOl) were "self-orthogonal" (as it has been accepted 
in [66]) , i.e. had zero binorms in ([9]), the would-be unity in (|10p would annihilate such 
states thereby signaling the incompleteness. 
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The PT-symmetry of a potential entails a related symmetry of eigenfunctions, 

V*(x) = V(-x) => Mx) = <(x) = 7 n^(-^), |7n| = 1- (11) 

when the PT-symmetry is not spontaneously broken, i.e. A* = A n (further on, for clarity, 
we restrict ourselves only with a case of the unbroken PT-symmetry although nearly all 
results can be generalized to the case of spontaneous PT-symmetry breaking with pairs 
of eigenstates having mutually complex conjugated eigenvalues). The normalization (jHJ) 
leads to the value 7„ = ±1 . One can see that the biorthogonality does not, in general, 
provide positive binorms of states related by the PT-symmetry, 

+00 

J ip*(-x)il> m (x) dx = 7„<5 nm = ±S nm , (12) 
—00 

bringing negative norm states in the PT-odd sector. 

For non-Hermitian Hamiltonians one can formulate the extended eigenvalue problem, 
searching not only for normalizable eigenfunctions but also for normalizable associated 
functions for discrete part of the energy spectrum. Some related problems have been 
known for a long time in mathematics of linear differential equations (see for instance, 
[68]) . 

Let us give the formal definition. 

Definition 1. The function ipn : i{x) is called a formal associated function of i-th 
order of the Hamiltonian h for a spectral value A n , if 

(h - A n ) i+ Vn, 4 = 0, (h- A„)Vn,i # 0, (13) 

where 'formal' emphasizes that a related function is not necessarily normalizable. 

In particular, the associated function of zero order ij) n o is a formal eigenfunction of h 
(a solution of the homogeneous Schrddinger equation, not necessarily normalizable) . 

Let us single out normalizable associated functions and the case when h maps them 
into normalizable functions H. 

Evidently this may occur only for non-Hermitian Hamiltonians. Then for any normal- 
izable associated functions i/} n ,i(%) an d ip n ',i'( x ) the transposition symmetry holds 

+00 +00 

J hip nti (x)ip n 'jt(x)dx= J ip nt i(x)hip n ^i>(x)dx. (14) 
—00 —00 

Furthermore one can prove the following relations: 

+00 

J t/j n ^(x)t/j n r ji f(x)dx = ,Vv,i') = 0, A n /A n /, (15) 
—00 

where (...,...) is scalar product. 
2 It takes place for a certain class of potentials described in Sec. 6, see Part II of our paper. 
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As well, let's take two normalizable associated functions tp n ,k(x) and ip n ,k'(x) so that, 
in general, k 7^ k' and there are two different sequences of associated functions for i < k 
and i' < k' 

ipn,i( x ) = ( h ~ A„) fc_ Vn,fc(a;), ipn,i'(x) = (h - Xnf ^n,k' {x) . (16) 

Then 

+00 

^n,iix)ij} ni ii{x)dx = ,Vv) = °> i + «' < max{fc, A:'} - 1. (17) 



In particular, for some normalizable associated function ijj n< i(x), the "self-orthogonality" 
661 is realized , 



ift, l {x)dx = 0, ^ n , l {x) = {h-\) i ~ l ^ n , i {x), 1 = 0,...,—— . (18) 



All the above relations are derived from the symmetry of a Hamiltonian under transposi- 
tion and the very definition of associate functions. 

We proceed to the special class of Hamiltonians for which the spectrum is discrete and 
there is a complete biorthogonal system {\ipn,a,i}, \ipn,a,i}} such that, 

h\lpn,afi) = K\ll>n,a,o), {h ~ K)\lpn,a,i) = \lpn,a,i-l) , (19) 

h ] \^n,a,p n . a -l) = K$n,a, Pn , a -l), ( h1 ~ KMn,a,p n , a -i-l) = \^n,a,p n , a -i) : 

where n = 0, 1, 2, ... is an index of an h eigenvalue X n , 

a = 1, . . . , d n is an index of a Jordan cell (block) for the given eigenvalue, A ra ; 
d n is a number of Jordan cells for A n ; 

i = 0, . . . , p nta — 1 is an index of associated function in the Jordan cell with indexes n, a 
and p n ^ a is a dimension of this Jordan cell. We have taken a general framework which is 
applicable also for matrix and/or multidimensional Hamiltonians. But the main results 
of this and the next sections are guaranteed only for scalar one-dimensional Hamiltonians 
with local potentials. 

We remark that the number d n is called a geometric multiplicity of the eigenvalue A n . 
For a scalar one-dimensional Schrodinger equation it cannot normally exceed 1 (but may 
reach 2 in specific cases of periodic potentials and of potentials unbounded from below). 
In turn, the sum ^2 a p n ,a is called an algebraic multiplicity of the eigenvalue X n . 

The completeness implies the biorthogonality relations 

(4>n,a,i\'>Pm,b,j) = ^nm^ab^ij , (20) 

and the resolution of identity 

+ 00 d n Pn,a — 1 

/ = EE E ivvvXiwi- (2i) 

71=0 (2=1 1=0 

The spectral decomposition for the Hamiltonian can be constructed as well, 

+00 d„ Pn,a — 1 Pn,a~ 2 

/i = EE[ An E \'&ri,a,i){i>n,a,i\+ E lV , n,o,t)(^n,a ) i+l| • (22) 
n=0 0=1 i=0 i=0 
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It represents the analog of the block-diagonal Jordan form for arbitrary non-Hermitian 
matrices |69j. 

If existing such biorthogonal systems are not unique. Indeed the relations (|19p remain 
invariant under the group of triangle transformations, 

Wn,a,i) = Yl a v\^n,a,j), 
0<j<i 

Wn,a,k)= E Pkl\^n,a,l), (23) 

k<l<p n ,a — 1 

where the matrix elements must obey the following equations, 

«ij = "i+l, 3+1 = a i-j, = a 00 0) 

Pkl = Pk+1, l+l = h-l+p n , a -l=P k -U Pn,a-1) Pp n , a ~l, Pn,a~l ¥= °. (24) 

The biorthogonality (|20p restricts the choice of pairs of matrices a and (3 in (|23p to be, 

j$ = 6T x . (25) 

This freedom in the redefinition of the biorthogonal basis is similar to Eq. ([S]) and it 
can be exploited to define the pairs of biorthogonal functions ip n ,a,i(x) = {x\^) n>a ,i) anc ^ 
ip nj a,i(x) = {x\tp n! a,i) in accordance with ([9]). However one has to take into account our 
enumeration of associated functions il> n ,a,i{x) vs. their conjugated ones ipn,a,i(x) as it is 
introduced in Eqs. (fl~9l) 

1pn,a,i( X ) = 4>n,a,p n , a -i-l( X ) = {^n,a,p n , a -i-l\x) . (26) 

Then the analog of Eq. ([9]) reads, 

+00 

/ ^x^^^dx = 8 Mj . (27) 

—00 

We stress that this kind of biorthogonal systems is determined uniquely up to an overall 
sign. 

In these terms it becomes clear that the relations (|17p have the meaning of orthogo- 
nality of some off-diagonal pairs in the biorthogonal system {|f/Vi,a,i)> |^n,oj)} as 

^n,a,i{x) = {h- A n ) Pn ' a_1_ Vn,a,p„, -l(«), 
i>n,a,j( x ) = ^n,a, Pn , a -l-j(x) = (h - X n ) J 1p n ,a,p n , a ~l (x) ■ (28) 

When comparing with specification of indices in Eq. (|17p one identifies 
Pn,a — 1 — J 1 h i *-* i' . In both cases k = k' = p n ^ a — 1 . Then the inequality (fT7|l 
singles out off-diagonal binorms, i < j — 1. From Eq. f|28|) it follows that in order to have 
all diagonal binorms non- vanishing it is sufficient to prove that at least one of them is not 
zero because 

+00 

J lpn,a,o{x)ll)n,a,p n , a -l{x) dx 
—00 

+00 

= J [(h- A n ) p "-"-V„,a, P „^i(a:) 

—00 
+00 

= J 1pn,a,i{ X )^n,a,p n , a -l-i(x) dx / 0. (29) 

— oo 



1pn,a,p n , a -l( x ) dx 
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The latter is necessary for the completeness of the basis because of the absence of biorthog- 
onal pairs of basis elements made of bound and associated functions in the diagonal res- 
olution of identity. If some of such pairs in resolution of identity (|2ip were biorthogonal 
then this operator would be at best a projector but not an identity. 

For a scalar one-dimensional Schrodinger equation the geometric multiplicity d n of 
the eigenvalue A n cannot normally exceed 1 . The latter possibility of non-degenerate 
eigenstates will be implied throughout this paper. Thereby in the rest of the paper the 
index a = 1 = d n will be omitted . 

It is certainly interesting to examine what are specific features of biorthogonal bases 
for PT-symmetric systems. Let us restrict ourselves by one-dimensional systems with 
nondegenerate spectrum of eigenstates and with the unbroken PT-symmetry, A* = A ra . 
First one can easily find that the functions ip ntPn -i-i(—x) are normalizable solutions of 
the initial Hamiltonian h and therefore can be decomposed into a linear combination of 
its basis, 

$n, Pn -i-l(-x) = tpn,i(.- X ) = ^ ln,ij4>n,j(x), (30) 

o<j<? 

where 7^ = J n ,i-j due to Eqs. conjugated to (fT9l) . By complex conjugation of Eq, (j30j) 
and its further convolution with 7 nj /_j (and after changing the sign x — > —x) one comes 
to the conditions on matrix elements, 

ln,l-ill,i-j . i = 0; 7n,o = !• ( 31 ) 
j<i<i J< 

The analysis of the biorthogonality relations (|20j) and (I27j) for the PT-symmetric basis 
(|30p leads to the conclusion that all numbers 7n,i-j are it real . Then one derives from 



(|31|) that 7„ i-j = and the elements of a Jordan cell basis, the eigen- and associated 

j<i 

functions are simultaneously PT-even or PT-odd depending on the sign of 7 ni o = ±1 , 

4>n,p n ~i-l{x) = 1pn,i( X ) = 7n,oVv(-z)- (32) 

Thus a biorthogonal basis of eigen-and associated functions exists which is fully compatible 
with PT-symmetry . As a consequence of (|32l) the integrals ip n fl(x)ip^ i pn _ 1 (—x) dx 
remain real. 

We remark that in general case the existence and the completeness of a biorthogonal 
system is not obvious (especially if the continuum spectrum is present [32]) and needs a 
careful examination. 



3 The origin of non-diagonalizable Hamiltonians 
is level confluence 

Let us demonstrate that one can generate a Jordan cell of a Hamiltonian in the process of 
coalescing levels of the Hamiltonian . In the simplest case one can consider a Hamiltonian 
depending on the parameter /i, 

K = -d 2 + v{x-v) 

with two eigenfunctions ipi^x; /i), 

h^ipi^ix; fj,) = Ai i2 Gu)^i i2 (x;^). 
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Assume that the levels Ai(/x) and A 2 (/x) coalesce for \i = 

^i(x; no) = ip2{x;iio) = <po(x), Xi(fJ-o) = A 2 (mo) = A . 
Let us also suppose that the functions q^(x; y) and ^jjj- (x; fj) are normalizable. Then, it 



is evident that 



[^-a 1|2 (m)]-^- 

(Vo - A o) -5— (x;yo) 



dV 
dji 



x'iM 



^1,2, 



(x;no) = [Ai(Mo) - A 2 Oo)]<A)(zO 



and, thus, the functions 



Or- 



Ai(Mo) - A 2 (^ ) 



y?i(x) 



(33) 



form a Jordan cell of the second order for the Hamiltonian ho 



h 



h tp = X <p , (h - \ )ipi = <p . 

Let us proceed now to the case with three coalescing levels of h^: 

[K ~ \j(fi)]i>j(x; fi) = 0, j = 1,2,3, 

ipi(x; no) = ip2{x;yo) = tp3(x;no) = <Po(x), Ai(/t ) = A 2 (// ) = A 3 (/u ) = A . 
Let us introduce the auxiliary functions 

T\\x;h) = [1 + x{n - HoMi(x;n), ipf\x; y) = ipj(x; fi), j = 2,3, 

(x;u) - — — (x:u) 

1,2, 



A'»-A^) 



where the constant yt is chosen so that the associated functions of the first order 
i>l\x;no) 



X[(fio) - A' 3 (/i ) 



+ 



Ai(Mo) - A' 3 (^ ) 



and 



A 2 (Mo) - A3O0) 

(cf. with (|33|) ) are identical!. When using these auxiliary functions one can obtain the 
canonical set of associated functions, 



<Po(x), ipi(x) = il>^(x;no) = ip^' (x; /j, ) 



<P2(x) 



2[X'M - A' 2 (// )] 

which form Jordan cell of the third order for the Hamiltonian /io = hn ■ 

h Q ipo = A0930, (ho - \o)<Pj = <pj-i, j = 1, 2. 



3 The constant n exists because in one-dimensional case the difference of normalizable associated functions of 
the first order for the same cigcnfunction is proportional to this eigenfunction.This freedom has been discussed 
in the previous section, see Eq. (|23f . 
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Thus, we have shown that the confluence of two (three) levels (of algebraic multiplicity 
1) leads to appearance of a Jordan cell of the second (third) order. The described con- 
struction is illustrated by an example in Subsec. 7.2. In the case of confluence of a larger 
number n of levels (of algebraic multiplicity 1) one can construct the canonical chain of 
n normalizable eigenfunction and associated functions in the same way. 

One can examine also the confluence of levels of different algebraic multiplicity. Let 
us restrict ourselves to the simplest case, when for the Hamiltonian ha level Ai(/i) of 
algebraic multiplicity 2 coalesces for fj, = /io with level X 2 (p-) of algebraic multiplicity 1: 

hf.ipw = Ai(/u)^i , [h^ - AiOu)]V>n = tpio, 

ip 10 (x;fio) = ^ 2 (x;no) = tp (x), Xi{fM)) = A 2 (mo) = A . 
Again we introduce the auxiliary functions 



i>i {x;n) = ipu (a;; /i), ip!p (x; fi) = 



Ai(^) - A' 2 (m) ' 

where one should choose h so that the associated functions of the first order 
ipf) (x; (jl ) = ipii(x;no) and 

.(1), , ^(x; Mo )-^(x^o) „ . . 

Ai(^o) - A 2 (/^o) AiC/io) - A 2 (/i ) 
are identical. Therefrom one can get the appropriate set of associated functions 

tpo(x), (pi(x) = ip^(x;no) = ijj^\x;fi ), 

MX) = 2[A 2 (, )-A' l( ,o)] ' (34) 
which form Jordan cell of the third order for the Hamiltonian ho = ha ■ 

h Q ipo = A0920, (ho - A )Vj = <Pj-i, 3 = 1, 2. 



4 Non-linear SUSY for complex potentials 

Supersymmetric Quantum Mechanics (SUSY QM) in one dimension represents a concise 
way for an almost isospectral transformation between two quantum systems |70| - [78] 
(see the reviews |79| - [84] ). Conventionally it can be built for a pair of Hamiltonians h + 
and h~ assembled into a Super-Hamiltonian, 

rr_(h+ \ _ / -d 2 + V 1 { X ) 

\ hr J \ -d 2 + v 2 (x) 

where the potential V(x) is, in general, complex. The (almost) isospectral connection 
between h + and h~ is realized by the intertwining relations , 

h + q + = q + h-, q-h + = h-q~, (36) 



J=-3 2 I + V(x), (35) 



10 



with q being components of the supercharges, 

«=(!! 5 o)- « 2 ^°- < 37 > 

The isospectral relations (j36j) result in the conservation of supercharges or the supersym- 
metry of the Super-Hamiltonian, 

[H, Q] = [H, Q] = 0, (38) 

In general, its algebraic closure is given, by a non-linear (deformed) SUSY algebra, 

{Q,Q}=V(H), (39) 

where V(H) is a function of the Super-Hamiltonian [3H [83], I85]. 

The relevant supercharges are supposed to be generated by N-th order differential 
operators with smooth coefficient functions w k (x): 

N 

q ± ^qN = Y. w k( x ) dk > w± = const =(Tl) N - (40) 

k=0 

We focus our analysis on the non-Hermitian Hamiltonians interrelated by complex super- 
charges, i.e. the supercharges with complex and smooth coefficient functions w^(x) . In 
this paper we choose connected to by means of * - transposition, qjj = (q~^Y ■ 
Evidently the Super-Hamiltonian of Schrodinger type is self-transposed (*-symmetric) as 
only the scalar potentials are under consideration. 

The algebraic structure of a Non-linear SUSY for local Hamiltonians is exhaustively 
determined by the following theorem, 

Theorem 1: on SUSY algebras with transposition symmetry 
Let us introduce two sets of N linearly independent functions 

4>ni x ) ( n = 1) " " " N) which represent complete sets of zero-modes of the supercharge com- 
ponents, 

q^t = 0, q N = (<?+)'. (41) 

Then: 

1) the Hamiltonians h^ have finite matrix representations when acting on the set of func- 
tions (fin( x )> 

h ± ^ = J2sL<t>Z, (42) 

771 

2) the SUSY algebra closure with Q = Q l takes the polynomial form, 

{Q,Q 1 } =det [EI- S + } E=H = det [EI- S-] E=H = V N (H). (43) 

The proof in [85] is based on the quasi-diagonalization of matrices S 1 * 1 , i.e. on their 
reduction to the Jordan canonical form S 1 * 1 which is block-diagonal. Such a diagonalization 
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can be realized by non-degenerate linear transformations f2 of the zero-mode sets which 
induce the similarity transformations of matrices S 1 * 1 , 

N N 

4>t = E h± 4>7 = E ^m, 

m=l m=l 

S ± = fi ± S ± (J2 ± )" 1 . (44) 

Evidently the canonical bases of zero-modes of intertwining operators form the set of 
(in general, formal - not necessarily normalizable) solutions and associated functions of 
the Hamiltonians h T . These elements of canonical bases of intertwining operator kernels 
are named as transformation functions. The proof in |85j can be easily generalized to the 
complex case with transposition symmetry being built up along the same scheme. 

For the polynomial SUSY algebra there is a possibility that the intertwining operators 
may be trivially reduced by a factor depending on the Hamiltonian, without any changes 
in the Hamiltonians h^, namely, 

q% = P{h ± )p± = P iP(hT), (45) 

where P{x) is assumed to be a polynomial and N > M + 2. Thus some of the roots 
of associated polynomials may not be involved in determination of the structure of the 
potentials. 

This problem of disentangling the nontrivial part of a supercharge and avoiding mul- 
tiple SUSY algebras generated by means of "dressing" can be systematically tackled with 
the help of the following theorem which can be also extended for complex potentials. 



"Strip-off" Theorem 2 . 

Let's assume the construction of the Theorem on SUSY algebras with transposition sym- 
metry. Then the requirement 

that the matrix S~ (or S + ) generated on the subspace of zero-modes of the operator q~^ 
(or q^j) contains m pairs (and no more) of Jordan cells with equal eigenvalues A; in each 
pair and the sizes Ski and k\ + 5ki (5ki being the size of a smallest cell in the l-th pair) 
is necessary and sufficient to ensure that the intertwining operator q^ (or q^) can be 
factorized: 

m 
1=1 

where p^ are intertwining operators of order M = N — 2 YlbLi = Sj=i % w hich can- 
not be decomposed further on in the product \^0) type. Herein kj for m + 1 < j < n are 
sizes of unpaired Jordan cells . 

Remark 1. The matrices S 1 * 1 cannot contain more than two Jordan cells with the same 
eigenvalue A because otherwise the operator A — h^ would have more than two linearly 
independent zero-mode solutions. 

Remark 2. This theorem together with the Theorem 1 entails the essential identity of the 
Jordan forms S~ and S + (up to transposition of certain Jordan cells). 

Remark 3. The supercharge components cannot be stripped-off if the polynomial Vn(%) 
does not have degenerate zeroes. The latter is sufficient to deal with SUSY charges non- 
trivially factorizable, but not necessary because degenerate zeroes may well arise in the 
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ladder (dressing chain) construction giving new pairs of isospectral potentials. 



Remark 4. In general, a Super-Hamiltonian may commute with several different super- 
charges. In this case few hidden-symmetry operators arise. The optimization of such a 
system of supercharges till one or two independent ones, their stripping-off and the mini- 
mal structure of a symmetry operator has been investigated in [85] in details for Hermitian 
Hamiltonians. In so far as in the quoted paper the transposition was used as a main conju- 
gation operation, all essential results of [85] remain to be valid also for complex potentials. 

The intertwining operators (supercharge components) can be formally factorized into 
the products of elementary Darboux operators. Let (f>J = cf>j, j = 1, . . . , N ba3 the basis, 
in which the S + -matrix (see Th. 1) has a canonical Jordan form and Xj is an eigenvalue 
of S + corresponding to the Jordan cell, to which <j>j belongs. Then adapting the Lemma 
1 from [85] for non-Hermitian Hamiltonians one can prove the following statements: 

1) for the supercharge components the factorization holds, in particular, 



In 



N 



1 • 



(47) 



where the Darboux operators 



r j =d + Xj{x), 



can be chosen so that 



r- ... r 1 (pN-j+l 



0. 



,N, 



j = l,...,N; 



(48) 



r N( r N) t + A ^ 



N- 
'- h N ; 



(49) 



2) the chain relations take place 

+ A m = r r( r D* + ^i = hi, 
(rf )*rf + Ax = h+ = h , 

3) the intermediate Hamiltonians hi, I = 1, . . . , N — 1 have Schrodinger form: 

hi = -d 2 + vi(x), vi(x) = xl+i(x) - x'i+i( x ) + A '+i = xh x ) + Xi{x) + k, 

Vi(x) = v (x) = xl(x) - Xi(z) + Ai> V 2 (x) = v N (x) = Xn( x ) +Xn( x ) + -V/v, ( 50 ) 
but, in general, with complex and/or singular potentials; 



4) the intertwining relations are valid: 

hi(r^ +i y = (r7, -,) f hi,+i , r7,-,h 



l = 0,...,N -1. 



(51) 



Let us introduce the generalized Crum determinants made of solutions of the initial 
Schrodinger equation for the Hamiltonian h + as well as of some of its formal associated 
functions, 



4>n(x) <f>' N ( X ) 
4>N-l(x) ^v-lO^) 



4-M 



<pN-j+l(x) <p' N _ j+1 (x) ... 4>%_^ +1 (x) 



(52) 



4 In what follows, in the notations for (f)J, we omit their relation to a Hamiltonian when it is only one of 
these sets which is used, in order to avoid too heavy indices. 
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Then in virtue of Eq. (|48p one finds the representation for the intertwining operators, 

4>n{x) 4>'n( x ) 



1 



r- ... v x 



Wj(x) 

and, consequently, 



<l>N-j+i(x) <f>' N _ j+1 (x) ... (ffyjy+^x) 
1 d ... & 



j = l,...,N (53) 



Wj+i(a) 

Wj(x) 



j = l,...,N-l. 



(54) 



Hence the intermediate superpotentials are uniquely determined by the chosen basis of 
solutions and formal associated functions of a given Hamiltonian h + for a given ordering 



, , [wj(x)/wj-x(x)]' w'Ax) w r j _ 1 (x) . 
XjS =~ ^ J ) " =-^tt+ J ' , J = 1,...,N, w (x) = l. (55) 

Thus, from Eq. (|50p one obtains the chain relations between intermediate potentials, 

Wj(x) 



Vj{x) - Vj-i(x) = 2x'j(x) 
and furthermore 



In 



Wj-l(x) 



1,...,N; 



Vj(x) - v (x) = -2[lnwj(x)]" , j = 0, 



N. 



(56) 



(57) 



It leads finally to the connection between the components of the potential in the Super- 
Hamiltonian, 

V 2 (x) = V 1 {x) - 2[lnw N (x)]". (58) 

The above set of relations is well-known from the Crum theory [86] • However here we 
have extended them including not only solutions of a Schrodinger equation but also its 
formal associated functions. 



Remark 5. Let's comment the sufficient conditions to preserve PT-symmetry under SUSY 
transformations, i.e. after intertwining with a Darboux-Crum operator . They consist 
in requirement for all transformation functions (all elements of the zero-mode basis of 
the intertwining operator) to be symmetric or antisymmetric in respect to PT-reflection . 
Indeed the PT-(anti)symmetry of basis elements entails the PT-(anti)symmetry of their 
Wronskian wherefrom one derives the PT-antisymmetry of W'/W and the PT-symmetry 
of V2 = V\ — 2(ln W)" . This is why under such conditions the above presented theorems 
and the forthcoming ones are certainly valid . 

5 Non-diagonalizable Hamiltonians and normal- 
izability of associated functions 

In this section we start examination of how the Darboux transformation may change 
the structure of a non-diagonalizable Hamiltonian with a Jordan cell spanned by a set of 
associated functions. We keep in mind that an intertwining operator may annihilate part 
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of them as zero modes. As a result we find the additive composition of Jordan cells for 
partner Hamiltonians mediated by a Jordan cell for the Hamiltonian mapping of the 
zero- mode subspace of the intertwining operator qj^. 

The rigorous results can be obtained with a specification of the class of potentials 
invariant under Darboux transformations. For such a class the asymptotic normalizability 
of associated functions at one of the infinities is preserved by Darboux transformations 
and will be described in certain lemmas and corollaries. When the normalizability on the 
whole axis is achieved we call the related associated functions as normalizable in general. 
The detailed mathematical proofs are given in part II, here we give only the general ideas 
of the construction as well as the formulations of the theorems and lemmas and the related 
corollaries . 

Let us investigate the Jordan structure of the Hamiltonians h . In what follows we 
restrict ourselves to the particular class of potentials: 

Definition 2. Let K be the set of all potentials V(x) such that: 

1) V(x)€C£; 

2) there are i?o > and e > (i?o and e depend on V(x)) such that for any |x| > i?o 
the inequality Re V(x) > e takes place; 

3) lmV(x)/KeV(x) = o(l), x — ► ±oo (this is sufficient to ensure the reality of the 
continuous spectrum); 

4) functions 

r— — \ 2 f \nx)l \V"(X)\ \ 

±Ro 

are bounded respectively for x > i?o and x < —Rq. 

Remark 6. The last condition is not very rigid: it is fulfilled (if x — ► +oo), for example, 
for potentials: 

1) V{x) =ax"<[l + o{l)\, V'{x) = a 7 xT- 1 [l + o(l)], 
V"(x) = a 7 ( 7 - l)x 7 " 2 + o(x 7 " 2 ), a > 0, 7 > 0; 

2) V(x) =V + ax-"' [l+o(l)], V'{x) = -a 7 a;^- 1 [l +o(l)], 
V"(x) = o 7 ( 7 + l)^ 7-2 [l + o(l)], F >0, a E C, Re 7 > 0; 

3) V(x) = ax a e bxl3 [l + o(l)], V'{x) = ab/3x a+0 - 1 e bxP [1 + o(l)], 
V"(x) = ab 2 p 2 x a+2l3 - 2 e bxl3 [l + o(l)], a > 0, b > 0, a G M, > 0; 

4) V(x) = V + ax a e~ hx0 [l + o(l)], V f (x) = -ab(3x a+ ^' 1 e - bxf, [l + o(l)], 
V"{x) = ab 2 p 2 x a+2f3 ~ 2 e- bxP [l + o(l)], F > 0, a e C, Re b > 0, a G C, /3 > 0. 

A similar statement is valid for x — > —00. Thus one can find in this class of potentials 
both representatives with purely discrete spectrum and the Hamiltonians with continuum 
spectrum . 

Remark7. This class of potentials only partially overlaps with PT-symmetric set of po- 
tentials investigated recently [9j-[17j,|23j, [23] . We call the related class of Hamiltonians 
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as softly non-Hermitian as they certainly don't involve complex coordinates in definition 
of asymptotic boundary conditions. 

The set K is closed under intertwining of Hamiltonians, that follows from the 

Lemma 1: on invariance of the potential set K . Let: 1) h + = —d 2 + V\(x), V\(x) £ 
K; 2) h~ = —d 2 + V%{x), V^z) G Cm.; 3) q~^h + = h~q^, where qj^ is differential oper- 
ator of Nth order with coefficients from C^; 4) each eigenvalue of S + -matrix of (see 
Th. 1) satisfies one of the conditions: either A < or Im A / 0. Then: 1) V^x) G K; 
2) coefficients of qjf belong to C|°; 5j h + q^ = q%h~ , where q^ = {q^ff , and moreover 
coefficients of belong to C-g also. 



Let us now analyze the normalizability properties of associated functions. 



Definition 3. A function f{x) is called normalizable at +oo (at — oo), if there is 
(R-) such that 



+oo R— 

f\mr*<«o (/i/wi^< + oo). (60 

R + — oo 

Otherwise f(x) is called non-normalizable at +oo (at — oo). 

Using the asymptotics of formal associated functions (see lemma 9, part II), one can 
show that: 

1) when the potential in h (by which, in what follows, we imply either h + or h~) 
belongs to the class K, one can show that any formal associated function of ra-th order 
of h, normalizable at +oo or — oo respectively, for the spectral value A, satisfying either 
A < or ImA / 0, can be decomposed as follows 

n 

^2 a m<Pm( x )> a 3,n = Const > «n,u/°> ( 61 ) 

j=0 

where ipj^^x), j > stand for either ipj^(x) or (pj^(x) and they form a sequence of 
associated functions normalizable at +oo or — oo respectively, 

h<Po,n = A<A),U> ( h ~ A )^j,U = Pj-i.U. 3 > 1 ( 62 ) 

Correspondingly, any associated function of n-th order of h, non-normalizable at the same 
+oo or — oo, for the same spectral value A can be presented as follows 

n 

J2 ( b m<P3,n( x ) + c i,n<Pj,n( x ))i ( 63 ) 

3=0 

where Cj,-\{ = const, either & n ,u 7^ or c^-j-j. 7^ and <pj ^ ^ {x) , j > form a sequence 
of non-normalizable at +00 or —00 respectively associated functions 

h$0M = Aft>,u> ( h ~ = &-i,U> J ^ ^ ( 64 ) 

2) for the Hamiltonian with a potential from the class K, there are no degenerate 
eigenvalues, satisfying either A < or ImA / 0, i.e. the eigenvalues, whose geometric 
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multiplicity exceeds 1 (the eigenvalues, for which there are more than one linearly indepen- 
dent eigenfunctions). Hence, for the Hamiltonian with a potential from K there is no more 
than one Jordan cell made of an eigenfunction and associated functions, normalizable on 
the whole axis, for any given eigenvalue A such that either A < or ImA / 0. 
Properties of associated functions under intertwining are described by the 

Lemma 2. Let: 1) the conditions of the lemma 1 take -place; 2) ip n {x), n = 0, . . . M 
be a sequence of associated functions of h + for the spectral value \: 

h + ip = \(f , {h + - \)ip n = ip n -i, n > 1, 

where either A < or ImA 7^ 0. Then: 

1) there is a number m such that < m < min{M + 1, N}, 

Qn^u = 0, n < m, 

and 

ipl = qNfm+l, I = 0,...,M -m 
is a sequence of associated functions of h~ for the spectral value A: 

h~ip = Xipo, {h~ - \)ipi = ipi-t, I > 1; 

2) if the function <p n (x), for a given < n < M , is normalizable at +co (on -co), 
then the function q^(p n is normalizable at +00 (on —00) as well. 

Corollary 1. Since h + is an intertwining operator for itself and both eigenvalues of 
its S + -matrix are zero, then if (p n (x) is normalizable at +00 (at —00), then tpj(x), j = 0, 
. . . n — 1 are normalizable at +00 (at —00) as well. 

Corollary 2. If there is an associated normalizable function <p n (x) of n-th order of 
the Hamiltonian h with a potential from K, for an eigenvalue A, which is either A < or 
ImA / 0, then for this eigenvalue there is an associated function ifj(x) of Hamiltonian h, 
normalizable on the whole axis, of any smaller order j: 

cpj = (h- A) n ~Vn, j = 0, ...,n-l. 

Corollary 3. Let </V/( x ) ^ e a canonical basis of zero-modes of the intertwining 
operator q^, i.e. such that S + -matrix (from the Theorems 1 and 2) has in this basis the 
canonical (Jordan) form: 

h+ <Pi,0 = X m,0> ( h+ ~ X i)<Pi,j = j = l,-..,ki-l, 

where ki is a rank of a Jordan cell for Aj. Then there are numbers kf^ and k^, < k^ ^ < ki 
related to the Hamiltonian h + such that for any i the functions 

v^-(^). i = o,...,/c+ a -i 

are normalizable at +00 or —00 respectively and the functions 

¥^0*0, j = kf u ,...,ki - 1 
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are non-normalizable at the same +00 or — 00 . Thus one can derive that the number 
of functions ip~j(x) normalizable on the whole axis is equivalent to min{fc^,fc^} and the 
number of functions ip~-{x) non-normalizable at both ends is given by ki — max{^, k^} . 

Independence of these numbers k^ ^ ona choice of the canonical basis, in the case of 
stripped-off , follows from 

Lemma 3. Let: 1) conditions of lemma 1 take place; 2) q~^ may not be stripped-off. 
Then any two formal associated functions of h + of the same order for the same spectral 
value X, when being zero-modes of qjj, are either simultaneously normalizable at +00 or 
simultaneously non-normalizable at +00. The same takes place at -co. 

We refer the reader for the proofs of the Lemmas 1 - 3 to the forthcoming Part II. 

6 Interrelation between Jordan cells in SUSY part- 
ners 

The first result on interrelation between Jordan structures of intertwined Hamiltonians 
and on the behavior of transformation functions at ±00 follows from the 

Lemma 4. Let us assume that: 1) conditions of lemma 1 take place; 2) {^j} o-nd 
Wtj} are canonical bases ofkeiq^ and ker q~j^ respectively; 3) q^ cannot be stripped-off; 
4) h is an algebraic multiplicity of the eigenvalue Aj of S + -matrix (see Ths. 1 and 2). 
Then for any i and j function ip^-{x) is normalizable (non-normalizable) at +00 if and 
only if iff k ._j_ 1 (x) is non-normalizable (normalizable) at +00. The same takes place at 
—00. 

Corollary 4. In order that for the level Aj the Hamiltonian h + does not have eigen- 
functions and associated functions normalizable on the whole axis and the Hamiltonian 
h~ has a Jordan cell of multiplicity i/_(Aj) spanned by eigenfunction and associated func- 
tions normalizable on whole axis (the same number z/_(Aj) measures the dimension of the 
subspace of non-normalizable zero- modes of q~) it is necessary and sufficient that 

ki = max{fe^,fe^} + v-(Xi) = \k~^ - k^\ + 

where kf^ {kfA are numbers of functions <pfj(x) normalizable at +cx) (—00) (see the 
previous Sec). Thus if there are no eigenfunctions and associated functions corresponding 
to the level lambda Aj of the initial Hamiltonian h + , normalizable on the whole axis and 
one wants to get the final Hamiltonian h~ with a Jordan cell of rank v-(\i) spanned 
by an eigenfunction and associated functions normalizable on the whole axis, one must 
choose transformation functions such that they contain v-(Xi) (and no more) associated 
functions of h + non-normalizable at both infinities . 

A more precise result on interrelation between Jordan structures of intertwined Hamil- 
tonians and on the behavior of transformation functions is given in the 
Index Theorem 3: on relation between Jordan structures of intertwined Hamiltonians 
Let us assume that: 1) conditions of lemma 4 take place; 2) v±(X) is the algebraic mul- 
tiplicity of an eigenvalue A of h^, i.e. the number of independent eigenfunctions and 
associated functions of h^ 1 normalizable on the whole axis ; 3) if A = \, where \ is an 
eigenvalue of S ± (see Th. 1), then n±(Xi) is a number of normalizable functions at both 
infinities among ipfAx), j = 0, . . . , fej — 1 and no(Aj) is a number of functions normalizable 
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only at one of infinities, among tfifAx), j = 0, . . . , ki — 1. Then the equality 

- n+(Ai) = i>-(Aj) - n_(Aj) 
ta/ces p/ace /or any i. Moreover ifrio(X) > /or some A = Xj, then for this Xj 

v+(Xj) - n+(\j) = u-(Xj) - n_(Aj) = 0; 
and if X is not an eigenvalue of S^ 1 but X < or ImX / 0, then z^+(A) = f-(A). 

7 An example of non-diagonalizable Hamiltoni- 
ans made by SUSY transformations 

7.1 SUSY system with Jordan cell of rank 2 

Let us start from the Darboux transformation of the free particle Hamiltonian (which is 
trivially PT-symmetric) , 

h + = -d 2 (65) 

and build an isospectral Hamiltonian (its SUSY partner) which is reflectionless |80, 85|[87] 
due to spectral equivalence to a free particle system, 

2 2 a{x - z)sh(2a(x - x )) - 2ch 2 (a(x - x )) 

h = —O — Iba r—— ; rr ; ttt^ , DDI 

[sh(2a(x-x )) + 2a(x-z)} 2 ' v ; 

a > 0, xq G R, z G C, Imz^O 
with the help of the intertwining operators 

qfh* = h± q f, q 2 =d 2 - -a 2 + q+ = (q 2 )\ 

W(x) = sh(2a(x -x )) + 2a(x - z). (67) 

If xo = Re z = the Hamiltonian h~ reveals PT-symmetry. Otherwise PT-symmetry is 
not realized although the energy spectrum remains real. 

The operator q 2 can be factorized into two intertwining operators of first order in 
derivatives, 

W'(x) 

Q2=%Qa, q a = d - ath(a{x - x )), q b = d - + ath (a(x - x )), (68) 

with the intermediate non-singular Hamiltonian of the ladder construction of Sec. 4, 

1 ch 2 (a(x — Xq)) ^ 

The canonical basis of q 2 consists of two non-normalizable functions: 

(x — z) 1 
tPo(x) = ch(a(x - x )), (p±(x) = — — sh(a(x - x )) + ^ch (a(x - x )), (70) 

h+tpo = X ifQ , (h + - X )tPi = <Po , X = -a 2 , S + = f ^° (71) 
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On the other hand the canonical basis of consists of two normalizable functions: 



(72) 



which form the Jordan cell for h corresponding to the level Ao: 
h~<p£ = Aq^o ' ( h ~ ~ X o)vf = ft > S ~ = 



A 
1 A 

In relation to factorization (|68p one can show that the zero-mode of becomes: 
sh(ct(x — xq)) (x — z) _ W(x) 



(73) 



^ a 2a 2ch (a(x — xq)) 4a ch (a(x — xq)) 

In turn the eigenfunctions of h~ for continuous spectrum read: 



a 1 , 

n— ( 74 ) 
2 <Po 



1 



kr 



fcel, h~tp(x;k) =k 2 ip(x;k). 



1 



1 + 



ik W'ix) 



W"(x) 



a 2 + k 2 W(x) 2(a 2 + k 2 ) W(x) 



Akx 



(75) 



One can check that eigenfunctions and associated functions of h obey the relations: 

+oo +oo 

>fc(x))" 4,-0, /rf(«)rf(x)*-l. -0, 



V>(x; k)ip(x; —k') dx = 5{k — k'), 



(76) 



where the last relation is understood, as usual, in the sense of distributions. 

One can also find that the functions <Pq(x), <pf{x) can be obtained by analytical 
continuation of ip(x; k) in k, 



lim 

k— »±ia 



2kdk 



lim [(A; 2 + a 2 )^(x; / 

k^±ia 



{{k 2 + a 2 )^(x-k)) 



• A /-e =ftX0 



, , . 1 =p 2az , 

( X ) ^0 W 



For this model resolution of identity made of eigenfunctions and associated functions 
of /i~ can be obtained by the conventional Green function method: 



+oo 



5{x-x')= I ^(x;k)ip(x';-k)dk + ip+ (x)<pf(x') + (pt(x)<p$ (x'), (77) 
or in the operator form, 



\^,k)(^,k\dk + IVoX^oI + 1-01)^1 



(78) 
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where the Dirac notations have been used: 

(x\ip,k)=il>(x;k), (x$,k)=i/>*(x;-k), (79) 

(a#o,i ) = (a#o,l) = (W^o( x )) ' ( 80 ) 

^- t |^i) = A |^i) J (/i-t - A )|^o) = h^\^,k) = k 2 \4,,k) (81) 

We stress that 1-00,1 ) are analogs of \ip n ,a,i) from Sec. 2. In this notations the biorthogonal 
relations take the form, 

(fotyk) = 8 jk , $,A#D,i) = <^i, o |^,fc) = 0, $,k\?P,k')=6(k-k'). (82) 
Accordingly, the spectral decomposition of ft, - can be easily derived, 

+00 

fc-= / fc 2 |V>, k\ dk — a 2 |V'o)('0o| — ct 2 \i>i) (4>i\ + |V ! o)(V'l|- (83) 



Let's remind that, for a given A = Aj, fcj is the number of zero-modes of q%, k^ , 
are the numbers of zero-modes normalizable at one of ±00 labeled by j, j, ^±(A) are the 
numbers of eigenfunctions and associated functions of normalizable on the whole axis, 
n±(A) are the numbers of eigenfunctions and associated functions among zero- modes tpf^ 
normalizable on the whole axis and no (A) is the number of eigenfunctions and associated 
functions among zero-modes <*pf- normalizable only at one end. They are defined in 
Sec. 5, 6 . For this model they take the following particular values (i = 0): 

ko = 2, fc+ u = 0, u + (X) = 0, k md = 2, i/_(A) ■- 



2, A — Ao, 
0, A^A , 



2, A — Ao, 
0, A/Aq, 



ra + (A) = 0; n o (A) = 0, n_(A) = 
and the Index Theorem holds, 

— n + (Xj) = U-(Xj) — n_(Aj) = 0. 

7.2 Coalescence of two levels 

The Hamiltonian h~ with Jordan cell for bound state (|66h is a particular limiting case of 
the Hamiltonian h~ with two non-degenerate bound states (of algebraic multiplicity 1). 
The former corresponds to the confluent case of the latter one. One again starts from the 
free particle Hamiltonian, 

h+ = -d 2 

and obtains its SUSY partner 

h~ = -d 2 - 16a 2 x 



2a/3 



sh (2a(x - x ))sh (2f3(x - z))-2ch 2 {a(x - x ))ch (2/3(x - z))+2sh 2 ((3(x - z)) 



[sh (2a(x - z )) + f sh (2/3(x - z))f 
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by intertwining with the operators 

W'(a 

W(x) " v " ' r ' ' 2 W(x 



where 



, W'(x) n . 2 ^ lW"(x) , , 

a 



W{x) = sh (2a{x - x )) + -sh {2ft{x - z)), 

P 

xq G E, Imz/O, a>0 (or - ia > 0), < /3 < 



21m z 



In the case ft 0, ft a the canonical basis of g 2 zero-modes (transformation functions) 
consists of 

¥±p{x) = ch.{k±p{x - £,±p)), k ±f3 = a±ft, £ ±/3 = 

so that 

fr + </>±/3 = A ±/ 399 ±/3 , A ±/3 = -/c 2 ^ = -(a ± /3) 2 . 

One can check that in the case ft ^ the function PF(x) is a Wronskian of (/? +| g(x) and 
ip-p(x) divided by ft and in the case ft = it is a product of —4a and of the Wronskian 
of ifQ (x) and y?^(x). 

For ft = both functions ip±p(x) coincide with </?q (x) from ([70]) . and as well \±p = Xq. 
In this case (p^ (x) is a linear combination of (p + p{x) and (f-p(x) in the following sense, 



<Pi- a 



1 

+ 4g^° 

/3=0 



g^(A +/3 - A_g) 

In the case /3 7^ the canonical basis of zero-modes of consists of 
One can prove that 



+00 

^M*-**^) (84> 



(this formula is valid also in the case ft = 0). The fact that the integrals (|84p vanish for 
/? = is in line with (|76p . It follows from (|84p that the normalized eigenfunctions of h~ 
have the form 



The eigenfunctions of h for continuous spectrum read, 

U2 . fl2 + 1.2 , _ 1^ 



kR, Jr^(x; fc ) = k 2 ip(x; k), (85) 
where the branch of */ (fc 2 + q 2 + /3 2 ) 2 — 4a 2 /3 2 can be defined by the condition 



V(fc 2 + a 2 + /? 2 ) 2 - 4a 2 /3 2 = k 2 + o{k 2 ), k -> 00 
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in the plane with cuts, linking branch points, situated in the upper (lower) half-plane. 
One can prove the following limits, 



lim [J{k 2 + a 2 + (3 2 ) 2 - 4a 2 (3 2 4>{x; k)} = J \ + - e =F(«"o+/?*)^ ( x ), 

k->±i(a+/3) y7T \j [3 a 



lim y(k 2 + a 2 + (3 2 ) 2 - 4a 2 /3 2 V(x; k)} = ^^£ A /4 " ~ e^ ax °- pz) ^-B(x). 
k— »±i(a— 0) yvr y p a 

In the limiting case j3 = the eigenfunction (p^ (x) and the associated function (ff(x) of 
h~ (see Eqs. ([72|) ) can be derived from ip±p(x) , 

¥>o ( s ) = ~ 2 iValim[^/P'ijj +l 3(x)} = 2y/a lim [y/j3ip„p(x)], 



(fi(x) = 2yfa lim , y 
Resolution of identity in the case (3^0 takes the form: 



^[VP{^(x)+ii> +f3 (x))] 



5{x - x') = ijj +f3 (x)ip + i 3 (x') + ip_[3(x)ip_p(x') + / ip(x;k)ip(x;-k)dk (86) 



and one can show that in the case a > 

l^tV^^+^ty) + ^_ /3 (x)i/>_ j 3(x / )] = (^(x^^x') + (p+(x)tp+(x') 

(cf. with ([771)). 

One can also check that the biorthogonal relations take place: 

+oo +oo +oo 

ip + /3(x)i/j-/3(x) dx = 0, y ip±p(x)tj}(x; /c) dx = 0, ^ il}{x;k)il){x;—k')dx = 5{k—k'). 

— OO — oo — oo 

(87) 

7.3 Symmetry operators 

For the Hamiltonian h~ of Subsec. 1 and 2 there exists the antisymmetric symmetry 
operator of the fifth order 

R 5 = q-dq+, R 5 h- = h~R 5 , R\ = -R 5 . 

In the case of PT-symmetry, xo = Rez = it anticommutes with PT-reflection 6pT, 
R 5 9 P T = -8 P TR 5 . 

The wave function (Pq(x), the associated function ip~^ (x) and the wave function of 
zero-energy bound state 



^(x;0)= ' 



| _ 1 W"{x) 
~ 2a 2 W{x) 
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are zero-modes of R$ because of ifQ^x) are zero-modes of q^ and q2^(x;0) = Const. In 
Subsec. 2 (in the case (3^0) the wave functions ip + p(x), ift-p(x) and the wave function 
of zero-energy bound state 



U 2 , o2 1 W"(xh 



2vr(a 2 - (3 2 ) 



are zero-modes of R$ because of ip±p(x) are zero-modes of q^ and q^ipix; 0) = Const. In 
both cases the spectrum of S-matrix of the symmetry operator (an analogue of S + from 
Th. 1) consists of and — (a ± (5) 2 . 

We notice that if k ^ then ip(x; k) and ij}{x\ —k) are linearly independent wave 
functions for the energy level E = k 2 . Simultaneously, tp(x;k) is an eigenfunction of the 
symmetry operator: 

R 5 ip(x; k) = ik(k 2 + a 2 ) 2 ip(x; k), 

or 

R 5 ip(x; k) = ik[(k 2 + a 2 + (3 2 ) 2 - Aa 2 (3 2 ]ip{x; k), 

Thus the zeroes of the eigenvalue of R§ are related to eigenvalues of S-matrix of R$ in 
accordance to |85|. 



8 Conclusions and perspectives: peculiarities of 
non-Hermitian Hamiltonians with continuous spec- 
trum 

In our paper we have investigated a bound state part of the spectrum for the class of 
potentials among which one can find also those ones with continuum spectrum. The choice 
of this class has allowed to keep all the analysis well below a possible continuum threshold. 
For this part of the spectrum the relationship between Jordan cells of SUSY partner 
Hamiltonians is firmly controlled by the Index Theorem 3 and it was well illuminated by 
an exactly solvable system with two coalescing bound states. We remind that the rigorous 
proofs of all new Lemmas and the Theorem 3 are postponed to the second part of this 
paper [53] . 

Meantime the approaching to the continuum threshold yields more subtle problems 
with normalizable eigen- and associated functions in continuum which may have zero 
binorm. As a consequence it may cause serious problems with the resolution of identity 
. This interesting problem we investigated in [32J. Here we would like only to draw 
attention to a class of models where the continuous spectrum is involved and elaborate 
the resolution of identity. 

Let us consider the model Hamiltonians 

h + = -d 2 , h~ = -d 2 + - — ^—77, Imz^O, (88) 

[x — z) z 

intertwined by the first-order operators q^: 

h±qf = qfh*, qt = Td- — . (89) 

x — z 

If Rez = these Hamiltonians are PT-symmetric and respectively their eigenfunctions 
possess definite PT- parities. 



24 



The eigenfunctions of h of continuous spectrum can be found in the form, 



ij)(x\ k) 



1 



1 



ik(x — z) 



Akx 



k £ M\{0}, h~il){x\ k) = k 2 ifj(x; k). (90) 



In addition, there is a normalizable eigenfunction of h on the lower end of continuous 
spectrum: 



t/j (X) 



(x - z) 

Isospectral relations between h + and h~ take the form, 

ikx -i p ^ 

iktp(x; k), k 7^ 0, 



2ir lim[ikip(x; k)], h V'o = 0. 

fc^O 



-1 



ipo(x), k = 0; 



q + i/j(x; k) = —ik 



Akx 



k^0, 



q +^ Q = 0, k = 0. 



The eigenfunctions of h satisfy the relations of biorthogonality, 

+oo +oo 



0. 



(91) 

(92) 
(93) 



(94) 



Resolution of identity made of eigenfunctions of h can be built as follows, 

5(x — x') = / ip(x;k)ip(x';—k)dk, 



(95) 



where the contour C must be a proper integration path in the complex k plane which 
allows to regularize the singularity in (I90p for k = circumventing it from up or from 
down. To reach an adequate definition of resolution of identity one can instead use the 
Newton-Leibnitz formula and rewrite (J95]) in the form 



*(*-*-)=( /+/>(*; WC;-*)* 



ip (x)ip (x') sme(x-x') 2 sin 2 (|(x - x')) 

1 ? 7V~ ^ 1 \F~i V' £>0 - ( 96 ) 

tts ir{x — x ) ire(x — z)[x — z) 

One can show |32j that the limit of the 3rd term of the right side of (|96|) (as a distribution) 
under e [ is zero for any test function from n L 2 (M) but the limit of the last term 
of the right side of ([96]) under e [ is zero only for test functions from n L 2 (]R; |x| 7 ), 
7 > 1. Thus for test functions from C|° n L 2 (M; |x| 7 ), 7 > 1 resolution of identity can be 
reduced to, 



6(x — x') = lim 



+00 



+ 



tp(x; k)ip(x'; —k) dk 



ipo(x)ipo(x') 



TTS 



(97) 
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and for test functions from H L 2 (M.) to, 

— e +00 

5(x-x') = hmj( / + / ^{x;k)^(x';-k)dk 



-00 e 



l-2sin 2 (|(x-x'))]Vo(x)Vo(^) }>■ (98) 



e|0 

1 

Decomposition (197j) seems to have a more natural form than (I98p . but its right side 
obviously cannot reproduce the normalizable eigenfunction 

Mx) ?c^nL 2 (R;\ x \-y), 7 >1 

because of the orthogonality relations (|94p . One can show that 

+00 

lim / —sin 2 (Ux - x'))^ 2 {x)ij (x') dx = lirn[e~ fea; Vo(a:')] = M x ')- (99) 
eio J ire 2 ej.0 



Thus it is the 3rd term in the right side of (|98p that provides the opportunity to reproduce 
ipo (x) and thereby to complete the resolution of identity. 
The spectral decomposition of h~ in this case reads 

k 2 \tp,k}(ip,k\dk, (100) 

—00 
where 

(x\ip,k) =ip(x;k), {x\$,k) =ip*(x;-k), h~^\^,k) = k 2 \^,k) . 

For the Hamiltonian h~ of this model there is an antisymmetric symmetry operator 
of the 3rd order 

3 3 
R 3 = q^dqf = -d 3 + -^d - , 

(x — z) A [x — z) 6 

R z hT = h'R 3 , Rl = -i? 3 , R 3 ip(x; k) = ik 3 ip(x; k). 

We notice that this model is a limiting case of the example of Subsec. 7.1 for a — > 
where z must be taken as a half sum of xq and z. 

One can generalize this model constructing the Hamiltonian (by intertwining with the 
Hamiltonian of a free particle) 

h- = -d 2 + ^l, 

[x — zy 

for which there is a Jordan cell, spanned by ffijr] normalizable eigenfunction and asso- 
ciated functions: 

n — 1 (2(n — 7) — 1)" 

^0 = 0, ^=^-1, ,- = (),...,[—], i> S {x)= mK2n _ mx _" z) n-2 j 

at the bottom of continuous spectrum . All these functions are mutually biorthogonal 
having also zero binorm. The problem of a correct resolution of identity seems to be 
solvable in the same way as for n = 1 . However the rigorous analysis is postponed to a 
forthcoming paper. 
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